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ABSTRACT 

The  engulfed  gas,  in  the  propagation  of  strong 
shock  through  a  y-law  gas  when  ^1,  is  confined  to  a 
thin  shell  near  the  shock  surface.   Based  upon  this 
picture,  several  models  are  constructed  for  the 
analysis  of  strong  shock  propagation  in  intiomogeneous 
media.   This  sequence  of  models  utilizes  successively 
more  detailed  pictures  of  the  wave  behind  the  shock. 
In  order  to  assess  the  validity  of  the  models,  com- 
parison is  made  with  a  numerical  perturbation  calcu- 
lation of  a  blast  in  an  atmosphere  with  linear 
variation  of  density.   It  is  found  that  the  distortion 
of  the  shock  front  and  the  distribution  of  mass  per 
unit  shock  surface  are  predicted  with  considerable 
accuracy  by  the  models,  whereas  the  transverse  mass 
flo^^^  is  sensitive  to  the  details  of  the  model  pictiore. 
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I.  Introduction 

A  very  large  air  blast.  In  the  several  megatons  of  TNT 
range,  spreads  throughout  such  large  volumes  of  the  atmosphere 
that  gravity  and  the  decrease  of  density  with  altitude  might 
be  able  to  Influence  the  blast  to  an  appreciable  degree. 
Exactly  what  this  influence  may  be  can  only  be  learned  by  a 
careful  numerical  integration  of  the  equations  of  motion. 
It  is  desirable,  nevertheless,  to  see  what  can  be  learned 
about  the  problem , by  several  more  qualitative  methods,  which 
give  somewhat  more  insight  into  the  problem. 

The  present  report  is  a  brief  presentation  of  a  few  of 
the  lines  of  approach  studied  by  the  authors. 

During  certain  times  in  the  motion  of  a  real  blast  it 
is  not  too  unlike  the  following  model,  which  we  study: 

(1)  The  blast  starts  at  a  point,  at  the  origin  of  coordinates, 

(2)  The  blast  is  "strong"  (in  the  sense  of  G.  I.  Taylor), 

(3)  The  atmospheric  density  (which  alone  is  of  consequence 
for  strong  shock)  decreases  linearly  with  altitude.   So  long 
as  the  shock  is  strong  the  hydro dynamical  accelerations  of 
the  blast  are  very  much  larger  than  the  acceleration  due  to 
gravity,  and  g  may  be  neglected  in  the  differential  equations 
of  motion.  Only  this  case  is  treated  in  the  present  report. 

Two  chief  lines  of  approach  are  discussed.   In  one  the 
hydrodynamical  equations  are  replaced  by  simpler  equations 
which  deal  only  with  integrated  flux  quantities,  according 
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to  the  method  of  the  excellent  "snowplow"  approximation. 
This  approach  gives  a  great  deal  of  insight  into  the  motion, 
and  very  likely  can  be  extended  into  the  region  of  weak  shock, 
Unfortunately  it  is  difficult  to  be  confident  of  the  correct- 
ness of  all  the  details  of  the  snowplow  approach  for  a  two 
dimensional  problem,  so  some  kind  of  standard  of  comparison 
is  provided  by  the  other  approach,  a  numerical  calculation  of 
the  blast  motion  by  the  method  of  perturbation  theory   • 
The  perturbation  calculation  is  described  first. 


(ij  See  also  "~ 

Hellund,  Emil,  Perturbation  Treatment  of  Blast  in 
Inhomogeneous  Media,  Argonne  National  Laboratory 
(Unpublished  Report). 
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II.  The  Perturbation  Calciilatlon 

Except  for  very  large  blast  radii  the  gravity  effect  may 
be  treated  as  a  small  perturbation.  An  analysis  on  this  basis 
Is  described  below,  the  unpertvirbed  solution  being  taken  to 
be  the  Taylor  result  for  strong  shock. 

As  we  have  noted,  so  long  as  the  shock  is  strong  the 
hydrodynamical  accelerations  of  the  blast  are  very  much 
larger  than  the  acceleration  of  gravity,  and  g  may  be 
neglected  in  the  differential  equations  of  motion.  Gravity 
thus  influences  the  blast,  for  strong  shock,  only  in  having 
produced  the  atmospheric  inhomogeneity  which  the  blast  must 
encounter.   In  terms  of  p, ,  the  density  Just  outside  the 
shock,  and  u  the  shock  velocity,  the  strong  shock  boundary 
conditions  give  for  density,  pressure,  and  velocity  just 
inside  the  shock  surface  the  results 

p  =  p^^  (y+i)/(y-i)  (1) 

p  =  2p^  uV(y+1)  (2) 

u^=  2tf/(Y+l)  (3) 

The  density  p,   which  enters  these  equations  is  approximately 

p-L=  PqiI  -  (rA)cos  9]  ik) 

where  the  scale  distance  X  is  about  8  kilometers,  this 
being  its  value  if  the  atmosphere  is  isothermal  and  at 
ordinary  temperature;  ©  is  the  angle  with  respect  to  the 


-  6  - 


vertical,   r  the  distance  from  blast  origin.   It  is  the 
second  term  of   (i;)  which  is  the  forcing  term  in  the  perturba- 
tion calciilation,  producing  all  the  departures  from  spherical 
symmetry.   It  may  be  shovm  that  the  gravity  effects  can  be 
relied  upon  to  be  small,  the  assimiption  of  the  present  calcula- 
tion, only  for  blast  radii  very  much  less  than  eight  kilometers. 

In  spherical  coordinates  the  differential  equations  of 
hydrodynamics  take  the  form 

«^^r  +  l<i,  [u^  +  u2]  -  1^  [  <r  (ru.)  -  <^  ]=  -ld£   (5) 
JF"   237   ^    **    r    57        ^  Pd^ 

<)uq  ^  u^  ^Uq  ^  u^  Hrxxg)^   -  1  d^ 
7t~   r  7^   r  ~7r       pr  3© 

(sin  ©  u^)] 


(6) 


r  sin  9  7©        ®         ^^^ 
(sin  ©  u^)].        (8) 


r  sin  ©  7©         ® 
These  are  now  to  be  solved  approximately,  subject  to  the 
conditions   (1)  -  (I4.),   and  to  conditions  at  r  =  0  which 

will  be  described  later.  For  the  spherically  symmetrical 

(2) 
problem,  about  which  we  perturb,  Taylor    found  the 

similarity  solution: 

(2)  Taylor,  G.  I.,  The  Formation  of  a  Blast  Wave  by  an 
Intense  Explosion,  RC-210  (Civil  Defense  Research 
Committee  of  Ministry  of  Home  Secxirity). 
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P  =  Pq  g(y)  (10) 

u  =  AR;^/2h(7)  (11) 

y  =  (r/R^)   ,  Rq  =  (|  At)2/^  (12) 

In  order  to  establish  a  correspondence  with  the  Taylor 
solution  we  now  introduce  for  the  present  problem  the 
similarity  variable 

y  =  (r/R)   ,  (13) 

where  R  is  the  radial  distance  to  the  shock  front,  as  a 
function  of  ©*  With  the  introduction  of  this  variable  the 
description  of  the  problem  has  been  transferred  from  the 
coordinates   (r,  ©,  t)  to  the  coordinates  (y,  ©,  t)» 

As  a  next  step  the  differential  equations  are  simplified 
to  the  case  of  first  order  perturbation  theory.   The  procedure 
is  straightfoirward  but  lengthy.  Each  hydrodynamical  quantity 
is  written  as  an  expansion  in  the  parameter  X  and  substi- 
tuted into  the  equations  of  motion.  The  zero  order  solution 
is  of  course  given  by  the  Taylor  blast  wave,  while  the  higher 
order  equations  are  all  linear  inhomogeneous •  Separation  of 
variables  can  be  achieved  in  the  first  order  part,  separation 
being  possible  if  the  perturbation  functions  are  given  the 

TTI  th 

form  P^(cos  9)  t™  f(y),  where  P^^  is  the  n    Legendre 
polynomial  and  f  a  function  to  be  determined.   In  the 
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present  discussion  a  detailed  study  will  be  given  only  for 

2/5 
the  solution  which  varies  according  to   t  '         and  according 

to  the  Legendre  function  of  first  order.   This  solution  will 
be  discussed  in  full.  -^ 

Written  out  explicitly,  oxir  solutions  for  pressure, 
density,  material  velocity  components,  and  shock  radius,  have 
the  following  form 

R 


p  =  y'^  Po  a^  r;3 


P  = 


u 


u   = 

9 


R  = 


A  R 


A  R 


'o 

-3/2 
o 

-3/2 

R 


f  +  -^  cos  ©  P(y)] 


g  +  -^  cos  ©  G(y)] 

R„ 
h  +  -~  cos  ©  H^(y)] 

R„ 

-^   sin  ©  H^Cy)] 

R 
1  +  _£  cos  ©  •  S] 


(li;) 
(15) 

(16) 
(17) 
(18) 


The  functions  P,  G,  H, ,  H^   are  determined  by  substitution 
in  (5)  -  (8)   and  now  satisfy  ordinary  differential  equations 

in  the  variable  y.   It  should  be  noted  in  the  above  that 

(3)      The  plausibility  of  the  restriction  of  the  study  to 
the  particular  solution  Just  mentioned  may  be  indicated 
in  the  following  way  J   Imagine  a  spherically  symmetric  a,l 
blast  moving  outwards  from  the  origin  at  time   t  =  0. 
Very  near  the  origin  the  second  term  of  equation   (I^.)  is 
negligible,  and  the  blast  remains  spherically  symmetric. 
The  spherical  symmetj^y  disappears  as  the  motion  proceeds, 
however,  as  terms  proportional  to  the  second  term  of   (ij.) 
gradually  appear  in  the  functions  which  describe  the 
blast.   That  these  perturbation  terms  are  proportional 
to  the  second  term  of  (i;)  follows  because  we  are  working 
with  linearized  equations.   Finally,  the  solution  described 
in  the  preceding  paragraph  is  precisely  that  one  which  has 
the  property  of  being  linear  in  the  forcing  term.  We 
sunnise  that  the  remaining  mathematically  admissible 
solutions  have  physically  inadmissible  singularities, 
but  shall  not  investigate  the  problem  at  this  time. 
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-r-  is  the  factor  which  varies  as  t  '  .  It  also  should  be 
noted  that  the  shock  distortion  is  measvired  by  the  constant 
S.  One  of  oxir  principal  tasks  below  will  be  the  determina- 
tion of  S, 

The  conditions  which  remain  to  be  satisfied  may  now  be 
listed.   First,  substitution  of  (ll|.)-(l8)  into   (D-d;) 
gives  the  boundary  conditions  at  the  shock  surface   (y  =  l)s 

H^(y=l)  =  U  S/(r+l)  (19) 

H2(y=l)  =  2  S/(y+1)  (20) 

F(y=i)  =  (|^)  ik   S-1]  (21) 

G(y=l)  =  -  (y+1)/(y-1)  (22) 

The  differential  equations  themselves  are  most  conveniently 
quoted  In  terms  of  the  set  of  variables 

^1  "^  %»   ^2  "^  ^2'   ^3  ^  -^/S*   ^k  ~   ^/S»        ^^^^ 
and  have  the  form 

dT  =  ^  ^ij  ^j  ^  ^^i    •  ^21,) 

The  coefficients  A^^   and  B,   are  complicated  ftinctions 
of  y,   involving  the  Taylor  blast  functions  f,   g,  h  and 
their  derivatives.  A  tabulation  of  the  coefficients  is  given 
in  the  Appendix.   In  the  present  work  a  solution  of  the 
equations   (2i|.)  was  obtained  by  the  very  straightforward 
procedure  of  step-by-step  numerical  integration,  the  integration 
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proceeding  inwards  from  y  =  1,  where  the  boundary  values 
(19)-(22)   are  known.  This  integration  was  carried  through 
for  the  case  y  ~  ^•k-t      and  ^o^  several  values  of  the  parameter 
S,  a  parameter  whose  meaning  must  now  be  explained* 

The  parameter  S  first  was  introduced  in  equation  (l8), 
and  entered  the  differential  equations  throu^  the  similarity 
transformation  (13)«  It  measures  the  distortion  of  shape  of 
the  shock  surface.  In  some  fashion  S  must  be  determined  by 
a  requirement  that  the  solutions  be  consistent  at  y  =  0» 
The  requirement  we  actually  use  is  that  S  be  chosen  so  that 
the  solution  is  "as  regular  as  possible"  at  y  =  0. 

Physically,  the  reasonableness  of  determining  the  shock 
shape  by  a  condition  at  the  center  may  be  understood  in  the 
following  way:   The  hydro dynamical  equations  are  local 
equations,  so  an  approximate  Integration  which  proceeds 
inwards  from  a  small  region  of  the  shock  surface  can  only 
determine  the  properties  of  the  voliome  behind  the  shock  with 
respect  to  the  properties  of  that  surface  region  from  which 
the  integration  has  started.  To  determine  the  overall  shock 
shape  and  velocity  the  various  different  surface  elements 
must  commxinicate  with  each  other,  by  way  of  signals  trans- 
mitted through  the  interior.  This  "communication"  is 
equivalent  mathematically  to  the  requirement  of  the  consistency 
with  each  other  of  the  solutions  which  emanate  from  the 
different  points  of  the  shock  stjrface.   In  the  present  problem, 
in  which  a  simplified  angular  behavior  of  the  solution  functions 
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Is  used  (equations  (Ii4.)-(l8) ),  the  consistency  requirement 
reduces  to  the  requirement  that  there  be  a  point  (or  a  region) 
near  the  origin  where  the  angle  dependent  terms  have  the  value 
zero*  With  only  one  adjustable  parameter  such  a  precise 
consistency  requirement  cannot  be  fulfilled  exactly,  and  we 
must  develop  a  procedure  for  satisfying  it  as  well  as  possible. 
In  particular,  since  the  density  becomes  extremely  small  in 
the  deep  interior,  wide  fluctuations  of  velocity  in  this 
region  will  be  physically  important  compared  to  density  and 
pressure  fluctuations. 

In  the  region  near  y  =  0   (a  region  which  extends  to 
y  values  as  large  as  0»8l )  the  equations   (Zlj.)  pass  over 
to  the  form 

5^  =  ^3^11^1+3  5^  .  (25) 

where         A^  >  =  lim  A.,  (26) 

■^J     J — >o  -^J 

B.  =  11m  B.  (27) 

1  y— >0  ^ 

These  simplified  coefficients  are  given  in  the  Appendix. 

It  is  very  easy  to  find  all  possible  linearly  independent 

solutions  of  equations   (2$).  Then  any  one  solution  of 

equations   (19)-(2i4.)  which  we  may  be  considering  must  go 

over  at  small  y  into  some  definite  linear  combination  of 

the  solutions  of  {2S),      and  its  behavior  near  y  =  0  may 

be  understood  in  terms  of  the  properties  of  these  elementary 

solutions. 
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The  elementary  solutions  of  (25)  have  the  form 

X,  =  0.70I4.  sy  +  2Z  ctn  y  ^  (28a) 

i 


1=1 

x^  =  O.Olii-  Sy  +  XZ  ^21  ^ 


(28b) 


.2   .  ^  _        ,^^h 


X3  =  1.783  Sy^  +  2Z  ^3^  7       ^  (28c) 

x^  =  7.57  S  +211  %1  ^"       ^  (28d) 


with 


p^  =  5.O8U 


p2  3  =  -0.61^5  +  1  1.92 


(29) 


P^^  =  -  9.0U 


Evidently  the  solutions  Involving  ^^,  p^,  (3.  all  become 
infinite  as  y  — >  0.  Strong  infinities  in  the  fvinctions 
X,  and  X.   are  especially  unacceptable ,  as  gx^  and  gx. 

are  the  pressure  and  density  perturbations.  On  the  other 

7  ^ 
hand,   g  varies  as  y   ,   so  only  the  solution  Involving 

pi   really  corresponds  to  physically  Impossible  pressure 

and  density  (so  energy)  distributions  within  the  blast. 

We  conclude  that  such  a  value  of  S  must  be  chosen  so  that 

as  the  corresponding  solutions  of  (2ij.)   are  carried  into 

the  region  y  ?::>  0  they  pass  over  to  the  forms   (28)  with 

In  detail,  the  requirement  that  the  Pi   solution  should 
not  occur  in  the  numerical  solution  is  met  by  forcing  a 
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certain  quantity  Q  to  vanish.  This  Q  is  the  scalar 
product  between  the  adjoint  of  the  characteristic  solution 
corresponding  to  p.  ,   and  the  difference  between  the  x^ 
and  the  particular  solution  of   (S5)» 

a  =  x^  +  0.331  ^2  +  1.539  -^   -  o.203(73c^)  -  1.915  yS  .  (30) 

We  find  Q  =  0  is  obtained  if  S  =  0,162.   It  should  be  noted 

that   (30)   (as  well  as  the  differential  equations  and  bovmdary 

conditions)  are  liaear  in  S,   so  that  the  correct  value  of 

S  may  be  interpolated  from  any  two  values  of  S  and  the 

corresponding  Q's.  As  a  nvunerical  check,  three  values  of 

S  were  employed  in  the  computation. 

Some  interpretation  of  the  result  S  =  0,162  is  possible 

by  brief  examination  of  equation   (18),   It  becomes 

R 
R  =  Rq[1  +  0.l62(-5^)  cos  ©]   ,  (31) 

or,  the  pertxarbation  of  the  shock  shape  is  only  16/6  of 
the  perturbation  of  the  ambient  density.  The  shock  remains 
more  nearly  spherical  and  at  constant  altitude  than  one 
might  at  first  suppose.   Presumably  this  result  is  a  conse- 
quence of  the  very  rapid  propagation  of  signals  across  the 
interior  of  a  Taylor  blast,  the  speed  of  sound  becoming 
strongly  infinite  as  y  — >  0,   so  that  the  different  surface 
elements  of  the  shock  very  easily  influence  each  other.   This 
same  observation  makes  it  less  surprising  that  our  approximate 
solution  gives  velocities,  x,   and  x^,  which  become  infinite 
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as   y  — >  0. 

The  results  of  the  perturbation  calculation  are  tabvilated 
in  Table  I,  and  plotted  in  Figures  1,  2,  3»  for  the  value 
S  =  .IS,   In  the  range  presented   (y  =  ,6  -  1.0),   these 
res\ilts  are  not  significantly  different  from  those  at  the 
more  accurate  value  of  S  =  .162« 
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y  |h(yJ 
(=U:Gr 

Hi(y) 
') 

H2(y) 

H^iy) 

=H2-Sh 

g(y) 
(=p:Gr 

G(y)   |f(y) 
?)           {=p:GIt: 

P(y) 

S 

1.00  .833 

.2667 

.1333 

.0000 

6.000 

-6.000 

1.166 

-.1+200 

.16 

.98 

.799 

.2730 

.091+8 

-.0330 

1+.066 

-3.1+35 

.955 

-.2732 

.16 

.96 

.766 

.2798 

.0596 

-.0628 

2.885 

-2.081+ 

.813 

-.1690 

.16 

.91+ 

.737 

.2866 

.0289 

-.0890 

2.120 

-1.365 

.711+ 

-.1395 

.16 

.92 

.710 

•2935 

.0012 

-.1123 

1.600 

-  .950 

.61+1+ 

-.1078 

.16 

.90 

.685 

.3001+ 

-.021+2 

-.1337 

1.232 

-  .695 

.592 

-.0666 

.16 

.88 

•  662 

.3071 

-.01+79 

-.1536 

.961+ 

-  .530 

.S^k 

-.0711+ 

.16 

.86 

.61+1 

.3137 

-.0701+ 

-.1729 

.762 

-  .1+16 

.525 

-.0600 

.16 

.81; 

.621 

.3199 

-.0921 

-.1911+ 

.608 

-   .31a 

.503 

-.0513 

.16 

.62 

•  602 

.3256 

-.1131+ 

-.2099 

.1+88 

-   .281+ 

.1+86 

-.01+1+3 

.16 

.60 

.581+ 

.3307 

-.1310+ 

-.2278 

.392 

-   .21+1 

.1+72 

-.0385 

.16 

.78 

.567 

.3352 

-.1551+ 

-.21+61 

.316 

-   .207 

.1+62 

-.0336 

.16 

.76 

.550 

.3389 

-.1767 

-.261+7 

.255 

-   .179 

.kSk 

-.0295 

.16 

.71+ 

.531+ 

.31+18 

-.1985 

-.2839 

.206 

-  .156 

.1+1+8 

-.0261 

.16 

.72 

.519 

.31+37 

-.2211 

-.301+1 

.165 

-  .136 

.kk2 

-.0230 

.16 

.70 

.503 

.31+1+5 

-.2/|)|8 

-.3252 

.132 

-  .118 

.1+38 

-.0203 

.16 

.68 

.1+88 

.310+0 

-.2700 

-.31+80 

.106 

-  .103 

.1+35 

-.0182 

.16 

.66 

.1+71+ 

.31+20 

-.2971 

-.3729 

.081+ 

-  .088 

.1+31+ 

-.0162 

.16 

.61+ 

.1+59 

.3382 

-.3271 

-.1+005 

.066 

-  .075 

.1+31 

-.01)|l| 

.16 

.62 

.1#+ 

.3323 

-.3607 

-.1+317 

.052 

-   .061+ 

.1+31 

-.0131 

.16 

.60 

.1+30 

.3238 

-.3993 

-.1+681 

.Oij.1 

-   .051+ 

.1+30 

-.0120 

,16 

Table  1:  Nxoinerical  Values 
of  Deviations  from  Taylor 
Plow  for  S  =  .16 


Fig.  I 

Amplitude  of  Velocity  Perturbation 
with  Taylor  Normalization 
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Amplitude  of  Density    Perturbation 
with  Taylor  Normalization 
(S=.I6) 


Fig.  3 

Amplitude  of  Pressure    Per)urbotion 
with  Taylor  Normalization 

(S=.I6) 


Ill,  The  Snowplow  Model 

A.  Preliminary  Results 

The  rather  ciunbersome  computations  of  the  perturbation 
method  yield  a  detailed  description  of  the  hydrodynamical 
system  under  consideration.   It  seems  reasonable  to  seek 
also  a  method  which,  while  yielding  less  detail,  may  be 
significantly  less  tedious  to  apply  and  more  physically 
direct;  in  brief,  to  calculate  for  a  model  of  the  system. 
Hints  as  to  such  models  may  be  obtained  both  from  otot 
detailed  numerical  results  and  from  general  physical  pictures, 
Thus,  it  has  been  observed  ^  in  numerical  computations  on 
inhomogeneous  air  blasts  that  the  shock  surface  pressure,  at 
a  given  instant  of  time,  tends  to  be  very  nearly  uniform 
under  a  wide  variety  of  conditions.   Since,  cmidely,  sxjrface 
pressure  gradients  produce  mass  flows  which  tend  to  cancel 
the  gradients,  such  a  result  is  not  strange.   If  we  take 
this  picture  literally,  then  the  pressure  correction  F  of 
(lij.)  will  vanish  at  the  boundary,  so  that,  according  to  (21), 
we  have  S=  »2^,      Thus  the  order  of  magnitude  of  the  surface 
distortion  is  correct,  but  such  a  pict\are  is  not  really 
adequate  for  the  strong  shock  region. 

We  now  consider  at  some  length  a  model  whose  justifi- 
cation is  not  hard  to  assess,  and  which  therefore  csui  be 
successively  corrected  with  corresponding  improvement 

C^) Unpublished  work  done  at  New  York  University,  1956  by 
H,  B,  Keller,  G,  L.  Lewis,  D,  L,  Ludwig, 
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in  the  results  obtained.   It  is  based  upon  the  following 
considerations.   In  a  strong  shock,  as  y  — >  1»  ^^^ 
available  compression  ratio  increases,  so  that  most  of  the 
gas  swept  up  into  the  shock  simply  gathers  at  the  shock 
surface,  and  travels  along  with  it.   Alternatively,  from 
(3)»  as  Y  — >  1»  the  siorface  gas  velocity  u  approaches 
the  shock  velocity  U.  The  not  effect  is  that  as  a  point 
blast  propagates  into  a  medium  for  which  y^l*  'the  engulfed 
gas  does  not  spread  into  the  interior  region;  it  is  merely 
"snowplowed"  into  a  thin  layer  behind  the  shock  surface. 
The  effective  removal  of  one  dimension  from  the  hydrodynamics 
introduces  obvious  simplifications  in  the  analysis  of  the 
system. 

Before  entering  into  the  analysis  proper,  let  us  ask 
what  we  expect  to  obtain  from  such  a  model;  clearly,  details 
of  the  wave  behind  the  blast  are  washed  out  in  this  picture. 
The  gross  properties  we  seek  will  then  be  the  surface 
distortion  factor  S   (from  which  all  s\irface  hydro  dynamical 
quantities  are  obtainable),  the  total  mass  M  per  \init 
surface,  and  as  one  aspect  of  the  hydrodynamical  structure 
behind  the  blast,  the  total  transverse  mass  flow  per  \anit 
STirface.  To  compare  with  the  perturbation  calculation, 
we  must  specify  the  interior  extension  of  a  unit  surface: 
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it  is  to  be  a  cone,  normal  to  the  shock  surface.   ' 

To  compute  the  pertinent  total  quantities  on  the  basis 

of  the  perturbation  solution,  we  observe  that  since  the 

shock  surface  is  a  displaced  sphere  given  by   (l8),   the 

center  line  of  the  cone  runs  from  the  displaced  sphere 

center  at  R  —  z  to  the  surface  point  R,   and  is  of 
o  A. 

length  R   (z  denotes  the  unit  vector  in  the  z  direction). 
Thus  if  Y  represents  relative  distance  from  the  center  of 
curvature,  we  have  for  any  P(r), 

^TOT  =  «o  /o  ^  ^^^  ^  <1-Y^  «o  I  ^>  ^Y       <^2> 
or  if 

P(YR)  =  F(Y)  +  I  cos  ©  5  P(Y)   ,  (33) 

P^Q^  =  R(l-S  I  cos  ©)/J  Y^(P(Y)  +1  cos  ©  5  P(Y))  dY 

2 
+  R^  cos  ©/J  Y^  (1-Y)  P'(Y)  dY       {3k) 

which  yields,  after  an  integration  by  parts, 

^TOT  "  ^^J^  ^  ^^^^   dY  +  2  cos  ©  /J  Y^  5  P(Y)  dY 

-  2S  I  cos  ©/J  Y(l-Y)  P(Y)  dY]   .  (35) 

(6)  Since  the  mass  density  rapidly  falls  to  zero,  the  fact 
that  the  evolute  is  not  generally  a  point  need  not 
concern  us;  in  the  special  distortion  we  consider,  the 
shock  surface  is  a  displaced  sphere  of  radius  R  ,   and 
so  the  evolute  actually  turns  out  to  be  a  point. 
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In  particular,  inserting  the  expression  (10)  for  p,  the 
surface  mass  density 

M(R,©)  =  P„  R  /i  [g(Y)  y^  -  2  I  S  cos  ©  (1-Y)  Y  g(Y) 


+  ?  cos  ©  G(Y)  Y^]  dY  ,  (36) 


R 


which,   after  ntamerical  integration,  results  in 

M  =  ^  Pjj  R   (1  -   .775  I  cos  ©)     .  (37) 

Finally,  inserting  the  expression  for  the  tangential  velocity 
vu  (R,)),   the  mean  transverse  velocity  becomes,  to  lowest 
order, 

u^  =  ^"^^^"^  =  ^^  sin  ©/J  g(Y)  Y^  [H2(Y)  -Sh  (Y)]  dY/ 

/Jl   g(Y)  Y^  dY  ,  (38) 

from  which  we  obtain 

0-1/2 
u^  =  -  .081^.  A  ii-Y-  sin  ©  .  (39) 

We  thus  find  that  in  the  analytical  expressions 

R 
R  =  Rq  [1  +  S  -j^  cos  ©]  (l|Oa) 

M  =  J  p^  R  [1-m  S  cos  ©]  (Ij.Ob) 


"t 


=  -  bA  R"^/^A  sin  )  (l+Oc) 


the  perturbation  calculation  gives  for  the  parameters  s  » 
m,  b 
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S  =  .ISZ,  m  =  .775,  b  =  .08I4.  (l|Od) 

In  the  succeeding  analysis,  several  simple  models  for  the 
hydrodynamics  will  be  set  up,  and  it  will  be  seen  in  what 
manner  the  values  these  models  give  compare  with  those  of 
(l+Od). 
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B»  The  Ideal  Snowplow 

We  first  analyze  the  case  in  which  y  — >  Ij   then  the 
mass  shell  has  zero  thickness  and  so  the  structure  behind 
the  blast  may  really  be  neglected.   In  order  to  determine 
the  parameters   S  and  m,  we  employ  the  integral  principles 
of  mass  conservation  and  conservation  of  longitudinal  momentum, 
together  with  the  fact  that  the  interior  pressure  is  a  constant, 
independent  of  9  (no  mass  is  present  to  maintain  a  pressure 

gradient).   Consider  then  a  unit  of  shock  surface,  containing 

-> 

mass  M;   as  it  propagates  with  velocity  U   it  absorbs 

mass   p,U  per  unit  time.   On  the  other  hand,  the  unit  surface 
expands  on  propagation,  the  particle  direction  being  that  of 
the  surface  normal.   If  the  surface  area  is  given  by  ^  , 
then  clearly 

V.O  =  [^(t4dt)  -  a(t)/<^(t)  U  dt,     (41) 

so  that  with  (X{t)   =1,  U   denoting  the  unit  vector  in 

-> 

direction  Tj   we  have 

d^/dt  =  U  ^•U  .  (1+2) 

The  mass  conservation  law  is  thus 

Since  we  have  agreed  to  allocate  all  spatial  quantities  to 
the  shock  surface,  and  thus  moving  with  it,  we  have  quite 
generally 
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P  =  d^'VF     .  (i4|.) 

Employing  {144),  ik3)     may  now  be  rewritten  as 

^.  (M  V)   =  p^    .  ikS) 

It  is  in  equation  ik^)     that  the  ambient  density  enters 
the  calculations,  leading  to  the  non-spherical  resvilts 
which  will  be  found. 

Next  we  utilize  the  longitudinal  component  of  the  law 

of  force  f  =  ma  for, a  unit  surface.  Both  the  surface 

2 
pressvire  p  =  p,U   and  internal  pressure,  say  p.,   are 

applied  longitudinally  i.e.  in  the  direction  of  the  surface 

A 

normal  U.   Furthermore  we  have 

U  ._d  tf*=  U  .  (u^.y )  U^  =  1  U^.  (U  .  V)  U^        (U6) 
dt  U 

=  1  (U^.V)  1  U^.  U*  =  U^.VU  , 
U         2 

whence 

M  u'^.VU  =  Pi  -  Pi  U^*  ikl) 

Alternatively,  eliminating  p  between  ikS)      and   (i|.7), 
there  results 

V.  M  u'^=  p^/U  (i^S  ) 

which  may  be  interpreted  as  the  force  law  for  the  material 
in  the  elementary  cone,  where  now  the  front  s\H»face  of  the 
elementary  cone  is  taken  to  be  on  the  ambiant  side  of  the 
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shock* 

Equations  ([|.5)  and  (1^.8)  are  stiff Iclent  to  determine 
m     and  S  In  (i|.0).  We  must  first  express  all  pertinent 
surface  quantities  in  terms  of  s.   Since 

^o  "  ^2  ^*^^^^  ^^2) 

R  =  R  (1  +  S  ^  cos  ©)  (18) 

O  A 

we  have,  retaining  only  first  order  terms  in  S, 

t  =  ^  R^/2  (1  -  I  s  S  cos  d)  (49) 

whence  u  =  Vt/|Vt|^  or 

t=   A[R"^/^  (1  +  I  S  I  cos  ©)r  +  R'-^/^  f  sin  ©  S]  (5la} 
U  =  AR"^/-  (1+1  S  I  cos  »)  (5lh) 

U=r+  s|sin©©  (5lc) 

at  the  sxirface  position  (R,©)»  Inserting  (^1)  as  well 
as  the  analytic  expressions  (i4.0a)»  (^Ob)  for  M  and  p, 
into  ikS)     and  (I1.8)  we  find 

p^  (1  -  I  cos  ©)  =  -^  (3-ipn  I  cos  ©)  +  I  p^  S  I  cos  ©  (52) 
^  r"'^  p^  (1  +  [|  S  +  ^  S  +  I  (|  S  -  m)]  5  cos  ©)  = 

=  Pi(t)  .  (53) 

But  p.,  a  function  of  t  alone,  must  then  be  oC  t  '-^,  or 
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2 
Pi  =  f  R"^  Pq  (1  -  3S  I  cos  ©)      t  =  I  Pq  ^^»   ^^^^ 

from  which  we  obtain 

l4Jn  -  2S  =  3  5in  -  23S  =  0  (55a) 

S  =  ^  =  .183  "  =  bI  "  '^^^   •  ^^^^^ 

Even  this  crude  picture  thus  yields  results  tolerably  close 
to  the  numerical  expression  (l+Od). 
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C,   Corrections  to  Longitudinal  Flow 

If  Y  is  greater  than  unity,  the  mass  shell  Is  no  longer 

of  zero  thickness,  and  so  some  aspects  of  the  structure  behind 

the  blast  make  their  presence  felt  even  in  otir  integral  type 

of  fonnulation#  These  corrections  are  necessary  even  for 

propagation  in  a  uniform  medium  (witness  the  exceptionally 

high  value  of  p^/p.lJ   in  (5U))  which  is  an  advantage, 

for  explicit  effects  of  inhomogenelty  may  be  omitted  in  their 

discussion*  Briefly,  the  thickness  of  the  mass  shell 

manifests  itself  in  that  the  effective  inert ial  mass  is  no 

longer  precisely  the  total  mass,  and  the  active  force  is  no 

longer  Just  the  difference  between  interior  and  surface 

pressure.  There  are  also  factors  to  be  inserted  in  the 
I 
relations  between  particle  and  shock  velocity,  and  surface 

pressiire  and  p-iTJ*^,  but  these  are  obvious,  and  indeed  their 

insertion  would  not  have  affected  the  results  of  the  previous 

section* 

Consider  first  the  mass  correction*   Our  statements  on 

mass  conservation  hold  vmchanged,  so  that  (i|.5)   remains 

correct  (at  least  for  the  homogeneous  case)*  Now  since 

the  mass  is  to  be  taken  as  moving  at  the  surface  particle 

velocity,  we  may  define  the  effective  longitudinal  inertial 

mass  Mj.  per  unit  surface  by 

p(r)  uj^  (?)  d^r  =  Mj^  ^«  V  u(R,©)   ,  (56) 


/- 


where  u-  is  the  longitudinal  component  of  the  material 
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time  derivative  of  the  particle  velocity,  and  the  integration 
is  over  the  cone  of  unit  shock  surface.  If  we  specialize  to 
the  homogeneous  Taylor  solution,  this  readily  yields 

a^^  S  ^  =  (y  +  1)/J  [|  h(Y)  +  (Y-h(Y)h'(Y)]  g{Y)y^   dY,  i$l] 

a  constant* 

Next,  the  longitudinal  force  on  the  cone  would  be 
p.  -  p(R,©)  if  the  pressure  were  to  drop  to  p^  at  once 
on  passing  inside  ,the  shock  surface*   Instead  the  quantity 
p  -  p.   is  rapidly  but  not  immediately  damped  out,  so  that 
presvimably  some  multiple  of  p--p(R,©)   should  be  subtracted 
from  p(R,©)-p^  Itself  to  produce  the  total  longitudinal 
force  Pj-i  we  write 

^L  ^  ^2  ^Pi  "  P^^»®^^  •  ^^^^ 

Since  the  total  force  is  obtained  by  a  sxirface  integration 
of  the  pressure  over  the  cone,  we  have 

O2  =  -  /p(i^)  U»d^{p^-p(R,©))   ,  (59) 

or  specializing  again  to   the  homogeneous  Taylor  case, 

ttg  =  [f(l)   -  2/J  Y  f(Y)   dY]/(f(l)   -  f(0))      ,  (60) 

again  a  constant. 

Now,  assxaming  the  effect  of  inhomogeneity  on  the 
values  of  a,   and  a^,  being  a  correction  on  a  correction, 
will  not  be  important,  our  conservation  and  force  equations 
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ik-S)     and  (1|.7)  become  Instead 


V  (M  U)  =  p  (61) 


2  .  „-?^.T-r  „  -  .   .-    2   _  .,2 


Y+: 


a^  M'^«V  ^  =  «^2  ^Pl  '  ^  Pi  ^    '      •  ^^^^ 


Solving  as  before,  we  find 


1*1 


Pi/pCRg)  =  1  -  I  ~     and  (63a) 

ipn  -  2S  =  3  (63b) 

3  ^m  +  (2U  -  5  ^)  S  =  6  (63c) 

2  2 

from  which 

2k   -  9(a^/a2)         _  84  -  l5(a^/a2) 
^  =  96  -I4(ai/a2)   '    °»  -  ^6  -  2hlo-2,/<i2^      *     ^^^^ 

If  we  adjust  the  ratio  a^/a^     to  correspond  to  the  Taylor 
value  of  •374  for  Pj^/p(Rq),  we  obtain 

a^/a^  =1.252       and  (65a) 

S  =  .163,  in  =  .831,  (65b) 

in  very  good  agreement  with  the  values  of  (l|.Od). 

It  must  be  emphasized  that  due  to  the  equations  of 
motion,  in  the  form 

p(?)u  j^(?)  +  Vl  p(^)  =  0  »  (66) 

relation  (62)  becomes  an  identity  on  insertion  of   (56) 
and  (59).  Thus  we  fail  to  be  surprised  when  on  numerical 
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computation  of     a,      and     a^     from     (57),      (60): 

a^  =  laoi,       a^  =  .879,     a^/a^  =  1.252,  (67) 

the  ratio  °-\/^2.     ^^  precisely  as  obtained  above.  The 
important  fact  ia  that  when  the  constants  a^, ,  a^     resulting 
from  this  identity  for  a  tmifonn  medium  are  carried  over  to 
the  non-imiform  case,  they  still  appear  to  carry  the  major 
part  of  the  thickness  connection  with  them. 
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D*  Other  Corrections 

We  have  seen  that  the  snowplow  model,  with  nominal 
corrections  to  the  longitudinal  flow,  yields  excellent 
agreement  with  the  shock  distortion  of  the  (presimiably  correct) 
perturbation  treatment,  as  well  as  a  good  estimate  for  the 
pattern  of  mass  flow.  Further  Improvement  may  be  expected 
as  the  Internal  structure  of  the  flow  Is  taken  further  Into 
account*  Thus,  since  the  Integrated  mass  xmlt  which  we  have 
been  considering  moves  on  a  curve,  there  must  be  specific 
transverse  effects  taking  place:  the  transverse  acceleration 
may  be  Imposed  by  transverse  pressxire  gradients  and  may  be 
reduced  by  a  transverse  mass  flow* 

The  transverse  pressure  gradient  Is  due  to  two  related 
effects.  First,  even  If  the  shock  surface  pressure  were 
constant,  the  radleJL  component  of  the  transverse  direction 
would  extract  a  non-vanishing  gradient  from  the  unpertxirbed 
Taylor  pattern.  Second  the  shock  surface  pressure  Is  not 
constant,  and  so  there  exists  a  transverse  gradient  which 
decays  on  passing  into  the  Interior.  The  mass  conservation 
equation  Is  also  modified  by  the  transverse  flow,  for  not 
only  is  there  a  direct  mass  flow  MU-,,  but  there  is  as 
well  an  Indirect  flow  due  to  the  change  of  the  elementary 
mass  cone  in  shape  and  orientation,  aside  from  that  of 
maintaining  its  unit  shock  surface,  >diich  has  already  been 
taken  into  account. 

If  one  Includes  the  aforementioned  transverse  effects 
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in  a  manner  similar  to  that  employed  for  the  longitudinal 
corrections  (gleaning  most  of  the  information  from  the 
original  Taylor  pattern)  only  the  mean  transverse  velocity 
parameter  b  is  required.  The  additional  transverse  force 
equation  wMch  is  now  available  then  enables  S,  m,   and 
b  to  be  determined*  The  value  of  b  thus  fo\md,  while 
always  small,  as  in  (i;Od},  is  quite  sensitive  to  the 
precise  formulation.   It  appears  then  that  a  somewhat 
haphazard  patching  up  of  the  original  model  to  accoxant  for 
these  additional  details  of  the  hydrodynamical  structure  is 
not  an  effective  way  to  proceed*  Rather,  one  should  set  up 
a  formalism  in  which  the  ideal  snowplow  approximation  appears 
at  an  early  stage,  but  which  may  then  be  systematically 
improved.  Such  an  approach  has  been  developed  and  will  be 
presented  at  a  later  date* 
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E»  A  Simplified  Approach 

It  Is  Interesting  to  point  out  that  one  can  describe  a 
very  simple  p3?oced\ire  for  obtaining  good  nvunerlcal  values 
for  the  parameters  S  and  m,   at  least  In  the  case  y  -   1»U» 
by  choosing  appropriate  aspects  of  the  snowplow  plctxire* 
Although  without  great  analytical  justification,  the  method 
is  readily  generalized  and  may  be  of  practical  value.  We 
reason  as  follows:  Neglecting  the  relatively  small  surface 
distortion  and  transverse  mass  flow,  the  variation  of  the 
mass  accumulation  at  a  unit  surface,  due  simply  to  the 

variation  of  ambient  density,  is 

2 
M(R,9)  =  /o  ^  Pi(r,©)  dr  ,  (68) 

or  inserting  (I|.), 

M  =  -| g—  i  cos  ©  ,  (69) 

whence  m  =  3/k-» 
If  we  now  append  the  simple  longitudinal  force  law   (I4.8) 
and  its  concomitant  in  (55a)  for  ovir  special  variation 
of  p, 

V*  M^  =  p^/u  ,        5m  -  23s  =  0  (70) 

we  obtain 

3  =  ^  =  .163.  (71) 
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P,  Concluding  Remarks 

Our  aim  has  been  to  Investigate  the  adequacy  of  the 
"snowplow"  type  of  model  for  analyzing  blast  propagation 
In  non-uniform  media.  As  a  standard  of  comparison,  the 
case  of  strong  shock  In  a  medivnn  of  y  =   l.l^.,   and  deviating 
only  slightly  from  homogeneity,  was  solved  numerically.   For 
this  special  case,  the  snowplow  with  longitudinal  corrections 
(Section  C)  and  the  empirically  constructed  model  (Section  E) 
reproduced  exceedingly  well  those  quantities  which  they  were 
suitable  for  determining:   shock  s\arface  position  (and  the 
surface  values  of  hydro dynamical  quantities  obtainable 
therefrom)  and  total  mass  per  unit  svirface,  Sinc«ji  neither 
the  strong  shock  boundary  conditions  nor  the  smallness  of 
the  Inhomogeneity  are  necessary  for  the  validity  of  the  model 
or  feasibility  of  the  computations  with  it,  it  appears  that 
such  models  may  be  valuable  tools  for  the  practical  solution 
of  shock  hydro dynamical  problems  for  intermediary  shock 
strength  and  for  large  density  and  pressure  inhomogenelties 
(with  Y  iio'*'  required  constant)  as  well* 
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Appendix 
Compilation  of  Coefficients  for  Equation  {2h) 

Define  D  =  f  -  gCy  -  h)^.  Then 


J;L  =  [^  i'  +  ^  ^'  +  g(7  -  h)  (|  h  +  2y  h'  )]  D" 


B2  =  f'/rsij'h) 

B3  =  [yf  (|  h  +  27  h' )  +  (y-h)  (3f+27f* )]  D"^ 
B^  =  [B^  -  2y(|-)]  (y  -  h)-^ 


A^2  =  -  2f/yD 

A^3  =  [2  -  Yli'  -  2Yh/y]  (g/rD) 


^U+  =  ^'<7  -  ^J/yD 


^21=° 

A22  =  (|  -  |)/(y-ii) 
A23  =  -  [Yy(y-h)]"^ 

^21,  =  ° 
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A32  =  -  2Yf    (y-h)/yD 
A3j^  =   (ff'/gD) 


A^  =  -  2g(y-h)/yD 

\3  "^  A^3(y-^)" 

Aj^  =   (fVvD)   +  (y-h)"^ 


In  the  limit   as     y  — >  0     the   above   go  over  to 

%  =  -  f   •        ^2  =  -  f  •        ^3  =  °   '        ^^11^  =  ° 

^31  =  ^  -  "^  •       ^32  =  2(l-r)  ,       A33  =  j^  .     131^  =  7[|  -  i) 
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